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ABSTRACT 

The FIRAS instrument on the COBE satellite has conducted an unbiased 
survey of the far-infrared emission from our Galaxy. The first results of this 
survey were reported by Wright et al. (1991). We report the results of new 
analyses of this spectral survey, which includes emission lines from 158 fim C + , 
122 /mi and 205 /mi N+, 370 /mi and 609 /mi C°, and CO J=2-l through 5-4. 
We report the morphological distribution along the galactic plane (b = 0°) 
of the spectral line emission, and the high galactic latitude intensities of the 
C + and 205 //m N + emission. In the galactic plane the 205 fim line of N + 
generally follows the 158 fim C + line distributions, but the intensities scale 
as J(N + 205 //in) oc J(C + 158 /mi) 1 ' 5 towards the inner Galaxy. The high 
galactic latitude intensity of the 158 fim fine structure transition from C + is 
I(C+ 158 /mi) = (1.43 ± 0.12) x 10" 6 esc |b| erg cm" 2 s" 1 sr" 1 for |b| > 15°, and 
it decreases more rapidly than the far infrared intensity with increasing galactic 
latitude. C + and neutral atomic hydrogen emission are closely correlated with 
a C + cooling rate of (2.65 ± 0.15) X 10~ 26 erg s -1 Hatom -1 . We conclude that 
this emission arises almost entirely from the Cold Neutral Medium. The high 
galactic latitude intensity of the 205 fim fine structure transition from N + is 
J(N + 205 //in) = (4 ± 1) x 10~ 8 esc \b\ erg cm -2 s -1 sr -1 arising entirely from 
the Warm Ionized Medium. We estimate the total ionizing photon rate in the 
Galaxy to be <f> = 3.5 X 10 53 ionizing photons per second, based on the 205 fim 
N + transition. 

Subject headings: ISM: atoms, bubbles, general, molecules — line: identification 
Galaxy: general 
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1. Introduction 

An unbiased far infrared survey of the spectral line emission from our Galaxy has been 
reported from the Cosmic Background Explorer (COBE) mission by Wright et al. (1991). 
Unlike other observations that concentrate on specific small regions of the sky at particular 
wavelengths, the COBE FIRAS instrument observed nearly the entire sky with nearly two 
decades of wavelength coverage. Thus, the COBE data show that the large scale diffuse 
emission from our Galaxy is dominated by a far infrared continuum that cools dust grains, 
and the 158 fim ground state transition of C + that cools the neutral gas in the Galaxy. 
Two transitions of N + trace the large-scale low-density extended ionized component of the 
Galaxy. In addition, the 370 fim and 609 fim lines of C° and the J=2-l through 5-4 lines of 
molecular CO were detected, and trace the interiors of neutral regions. 

Mechanisms for heating interstellar gas include collisions, radiation from stars, shocks, 
and cosmic rays. Examination of the spectral lines that cool the gas can help determine 
the dominant excitation mechanisms and conditions. Studies of particular bright regions 
in our Galaxy and observations of external galaxies have suggested that stellar ultraviolet 
radiation can ionize vast volumes of a galaxy and that the far ultraviolet (FUV) radiation 
impinging on neutral cloud surfaces is responsible for a large fraction of the observed far 
infrared (FIR) spectral line emission that cools the gas. 

Tielens & Hollenbach (1985a) defined photodissociation regions (hereafter referred to 
as PDRs) as neutral "regions where FUV radiation dominates the heating and/or some 
important aspect of the chemistry. Thus photodissociation regions include most of the 
atomic gas in a galaxy, both in diffuse clouds and in the denser regions..." Early work 
on PDR physics includes Glassgold & Langer (1974, 1975, 1976), Langer (1976), Black 
& Dalgarno (1977), Clavel, Viala & Bel (1978), and de Jong, Dalgarno & Boland (1980). 
Tielens & Hollenbach (1985a) present a detailed and quantitative model for one dimensional 
PDRs. 

The basic photodissociation process can be summarized as follows. Far ultraviolet 
photons in the energy range 6-13.6 eV, often from and B stars, impinge on the surfaces 
of neutral clouds where they will be stopped by impacts with dust grains. The dust 
grains emit photoelectrons and cool by continuum far infrared emission. The UV generally 
photodissociates molecules and photoionizes atoms with ionization potentials less than the 
Lyman limit. The photoelectrons ejected from dust grains heat the gas, resulting in the 
excitation of C + atoms from their ground state to their first excited state, AE/k = 91 
K above the ground state. C + cools by the emission of 158 fim radiation. The entire 
process thus converts far UV radiation to far infrared continuum and spectral line emission. 
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Only far UV radiation can contribute to this process since the typical work function of 
an interstellar grain is ~6 eV (de Jong 1980). The photoelectric mechanism can typically 
convert a maximum of ~4 % of the FUV energy into gas heating. This includes the 
assumption of a ~ 10% probability that an absorbed photon results in an ejected electron 
(the photoelectric yield). Wright et al. (1991) determined, based on COBE data, that the 
158 fim C + line alone emits 0.3% of the bolometric far infrared (ss UV) luminosity of our 
Galaxy. 

The UV flux incident on a neutral cloud, Go, is usually expressed in units of the Habing 
(1968) flux of 1.6 X 10~ 3 ergs cm -2 s -1 . For unusually highly excited regions, such as Orion, 
the UV flux incident on the neutral PDR region is Go ~ 10 4-5 and the gas density is n ~ 10 5 
cm -3 . The ionization potential of carbon is 11.3 eV, so carbon can be singly ionized where 
hydrogen is neutral. The ionization potential of oxygen is 13.6 eV, like that of hydrogen, 
so oxygen is neutral where hydrogen is neutral. In general, for high UV fields (Go > 10 2 ' 5 ) 
and/or high densities (n > 10 3 ' 5 cm -3 ), the 0° fine structure lines at 63 fim and 145 fim 
contribute significantly to the cooling of the neutral gas. For lower densities and UV field 
strengths, the 157.741 /mi (Cooksy, Blake, & Saykally 1986) transition of C + ( 2 Ps/2 — 2 -P1/2) 
is the major cooling line. This is because in 0° the 3 Pi excited state is 228 K above the 3 P 2 
ground state, and the 3 P state is 326 K above the ground state. In C + the 2 P3/2 excited 
state is only 91 K above the ground state. Dust grains become positively charged, especially 
near the surfaces of clouds, resulting in a decrease of the photoelectric heating efficiency. 
An increased electron density results in a higher dust-electron recombination rate, allowing 
the UV energy to be converted to gas heating more efficiently. Thus, the influence of an 
increase in density is similar to an increase in the gas excitation. 

The visual extinction into a neutral cloud, Ay } is related to the hydrogen column 
density, N H , by N H /A V = 1.8 x 10 21 cm" 2 mag" 1 (cf. Burstein, & Heiles 1978, Bohlin, 
Savage, & Drake 1978), although this is somewhat dependent on the local gas-to-dust ratio. 
Hydrogen in the PDR is predominantly in the form of H I for Ay < 2 and H 2 for Ay > 3. 
C + dominates the cooling of the outer regions of the cloud (Ay < 4) for low-excitation 
PDRs, as described by Hollenbach, Takahashi, & Tielens (1991). C° should dominate the 
cooling for 4 < Ay < 6 and CO rotational transitions dominate the cooling for Ay > 6. 
The Ay values, above, vary according to the detailed physical properties of the particular 
PDR (i.e. n, Go, and geometrical effects). Significant UV shielding (i.e. Ay ~ 10 — 20) is 
required before OH, H 2 0, and 2 can be sustained. For Ay > 6 heating from cosmic rays 
and gas-grain collisions become significant relative to the UV heating. 

The photodissociation model has been successfully applied to Orion, which is a face-on 
PDR (Tielens & Hollenbach 1985b), and to M17SW, which is an edge-on PDR (Meixner 
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1992). Wolfire, Tielens, & Hollenbach (1990) applied the model to the nuclear region of our 
Galaxy and to M82 to derive UV field strengths and gas densities. Schilke et al. (1993) 
also examined M82, but they find that the neutral atomic carbon emission is anomalously 
enhanced above the predicted levels from the standard PDR model. This may be due to 
an unusually high cosmic ray rate in M82. Low density PDR models, which are important 
for interpreting the average PDR conditions in our Galaxy, were computed by Hollenbach, 
Takahashi, & Tielens (1991). Wolfire, Hollenbach, & Tielens (1992) conclude that CO 1-0 
emission from giant molecular clouds is primarily the result of the PDR physics. Crawford 
et al. (1985) and others have found that PDR models account for the excellent correlation 
both between 158 fim C + emission and CO 1-0 emission, and between C + and far IR 
continuum emission. 

2. Observations &c Previous Results 

Results are presented from the COBE Far InfraRed Absolute Spectrophotometer 
(FIRAS) instrument. The instrument is a scanning Michelson interferometer that covers the 
wavelength range from 100 fim to 1 cm. While the primary purpose of the instrument is to 
make precise measurements of the cosmic microwave background radiation, the instrument 
design offers many advantages for the study of large scale diffuse emission from our Galaxy. 
For example, the FIRAS has a 7° beam, has the ability to make absolute measurements 
of intensity with no beam switching, and has observed nearly the entire sky. For these 
reasons, the FIRAS is well-suited for measurements of diffuse emission. The FIRAS spectral 
resolution is low, as discussed below in §3, and no spectral lines are resolved. Despite this, 
the general galactic rotation can be discerned using the detected FIRAS spectral lines. 

The first cosmological results from the COBE FIRAS instrument were reported by 
Mather et al. (1990), and recent overall summaries of COBE results are found in Boggess 
et al. (1992) and Bennett et al. (1992a, b). More recent cosmological results are reported 
by Mather et al. (1994), Fixsen et al. (1994a), and Wright et al. (1994). 

The first results of the COBE FIRAS unbiased far infrared survey of the Milky Way 
were published by Wright et al. (1991). Table 1 in that paper summarizes spectral line 
detections from COBE FIRAS: 158 //m C+, 370 //m and 609 //m C°, J=2-l through J=5-4 
CO, and 122 /mi and 205 fim N + . Since that time, improved data processing algorithms 
and techniques have allowed for a more detailed examination of the data. In addition, since 
the time of the Wright et al. (1991) analysis, a laboratory detection has been reported 
for the 205 //m ground state transition in N + . The frequency of the N + transitions are 
measured to be 1461.1319(6) GHz and 2459.3801(4) GHz (Brown & Evenson 1993). 
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The C + , C°, and CO emission lines detected by COBE can immediately be recognized 
as the expected cooling lines from PDRs with low densities and small UV fields (Hollenbach, 
Takahashi, & Tielens 1991). The result that the galactic-averaged PDR properties are 
described by lower excitation conditions than objects such as Orion is reinforced by the 
weakness of the 145 fim 0° line. 

The two transitions of N + are recognized as arising from the ionized medium. As 
pointed out by Wright et al., this emission can not be coming only from classical, 
high-density, H II regions. The ratio of the intensities of the 3 P 2 to 3 Pi 122 /mi N + 
line to the 3 Pi to 3 P 205 fim ground state transition is a probe of the density of the 
ionized gas. Rubin (1985) provides a grid of classical H II region models that result in 
J(N + 122 /im)//(N + 205 fim) ranging from 3 (for n e 10 2 cm -3 ) to 10 (for n e 10 5 cm -3 ). 
In the low-density limit, n e <C 100 cm -3 , the expected intensity ratio is 0.7. Wright et al. 
reported that the measured COBE ratio is 1.0 to 1.6, depending on the adopted values 
for the instrumental resolution at the two wavelengths. Thus the ionized gas measured by 
COBE FIRAS is likely to arise mostly from large, diffuse, regions. Reynolds (1990, and 
references therein) has extensively studied the warm ionized medium (WIM), concluding 
that a large fraction of the Galaxy contains ~ 8000 K fully ionized gas with a density of 
n e 0.2 cm -3 . However, the galactic disk is opaque to the Ha and other optical lines used 
to study this medium so only a local (typically < 2 kpc) view exists. The galactic disk 
is transparent to microwave continuum free-free emission, but it is difficult to separate 
the dust and synchrotron components from the overall microwave emission (Bennett et al. 
1992c). Radio recombination line observations have not been sensitive enough to observe 
the WIM. 

3. Line Fitting and Calibration 

The FIRAS instrument and its calibration were described in detail by Fixsen et al. 
(1994b). In this paper we use data from all four detectors (Left Low, Left High, Right Low, 
and Right High) and all of the interferometer mirror scan modes used in flight (Short Slow, 
Short Fast, and Long Fast). Data from all of the channels, scan modes, and frequency 
bands were combined using a weighted average to form a single data set with 167 spectral 
points from 4.4 mm to 104 /mi (2.3 to 96.3 cm -1 ). Right High Short Slow (RHSS) mode 
data dominate the high frequency spectrum (20-96 cm -1 ) and Left Low Short Slow (LLSS) 
mode data dominate the low frequency spectrum (2-20 cm -1 ). 
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The ideal instrument line response function for an interferometer is 

//'CO 
dO, J dv cos(27rz/x cos 6 + </>), (1) 

where x is the path length difference in cm between the two arms of the interferometer, v is 
the line frequency in waves per centimeter, 6 is the angle of each ray from the beam axis, 
dft is a differential of the solid angle of the beam, <f> is the phase shift of the interferometer 
at the frequency v } and Qbeam normalizes the solid angle integral over the beam. The 
instrument produces interferograms that we Fourier transform to recover the approximate 
sky spectrum, 

dxe 27Ttl/x A(x)I(x) (2) 

- CO 

or 

/CO P fGO 

dxe 27Tll/x A(x) dtt dvco^irvxcosO + <j>). (3) 
-oo J J0 

The factor A(x), the apodization function, is a Fourier window function that accounts for 
the fact that data are not obtained with equal weight over the full range of the integral. 

With A(x) = 1, the integral over solid angle (with cylindrical symmetry) produces a 
rectangular spectral line profile of width Av / v = 1 — cos(# ma:E ), where 9 max = 7° is the 
cut-off angle of the FIRAS interferometer. The effect of this cut-off angle can be represented 
in terms of an additional Fourier window, or self-apodization function, A se if(x), while taking 
cos(6 max ) = 1. Ideally A se ij(x) would be the Fourier transform of a rectangular spectral line 
profile, i.e. a sine function, but in practice there is some taper in the instrument profile that 
makes a Gaussian a better approximation than a sine function. For the FIRAS instrument 
Avjv = 0.007 FWHM by design, but the measured A se ij is approximately Gaussian with 
a FWHM of Avjv = 0.0035 at the C + line. A possible explanation for the discrepancy is 
that the expected large optical aberrations reduce the optical efficiency near the edge of the 
beam. 

For the FIRAS instrument, A(x) = where data are not taken, so the instrument 
line profile width is further broadened. Combining this effect with the self-apodization 
function, we take the total A tot (x) = A(x)A se ij. We take A(x) to be a particular asymmetric 
smooth window function to enhance the signal-to-noise ratio for wide spectral features and 
minimize ringing from spectral lines. The Fourier transform of A(x) is the dominant factor 
in determining the instrument line profile. The piecewise expression for the apodization 
function A(x{) for the 512 sampled data points is 

2 

fori G [1,512] (4) 



Ai — fi 
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where f 1 =f 2 = 0. For Left and Right Low Long Fast (LLLF and RLLF) data j 
c = 88, and 



513, 



fi 







1 



cos 



COS 



7r(i-2) 
30 



7r(i'+32-2c) 
30 



For all other scan mode data j ' = 1 and 
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" (2C ~3o 82 ~ 8) )}/ 2 



7r(513-i) 
30 



(6) 



(as seen in Figure 7 of Fixsen et al. 1994b), where for the Right and Left High Long Fast 
(RHLF and LHLF) modes c = 353, for all Short Slow scan modes c = 355, and for Right 
and Left High Short Fast (RHSF and LHSF) modes c = 356. 

Since A(x) is asymmetric, the resulting line profile is complex. The real (complex) 
part of the transform of A(x) has a width of 0.65 cm -1 (0.52 cm -1 ) for Short mirror 
scan mode data and 0.16 (0.13 cm -1 (0.13 cm -1 ) FWHM for Long mode data (usable 
only for frequencies below 20 cm -1 ). Thus the instrument spectral profile model is 
generated from the ideal cosine wave interferogram with the two window functions applied, 
Atot(z) = A se if(x)A(x). The resulting spectrum is calibrated, which involves dividing it by 
a frequency-dependent gain function. The final line profile is asymmetric and complex. 

Continuum and line spectra were fitted to the data as follows. A 2.726 K Planck 
spectrum was subtracted from the spectrum at each sky pixel. Complex parameters of the 
fit (and their corresponding degrees of freedom at each pixel) were: the amplitudes at each 
assumed line frequency (32), the amplitude of a d(Planck)/dT term to remove the cosmic 
microwave dipole (2), and coefficients of Legendre polynomials for the continuum emission 
of the dust (40). The fit removes 74 degrees of freedom, leaving 260 degrees of freedom at 
each of 5856 pixels. The data were weighted by an inverse variance that is dominated by 
detector noise, but includes other error terms (see Fixsen et al. 1994b). The use of this 
diagonal matrix weighting is an approximation to the full covariance matrix, which is not 
perfectly diagonal because of the effects of the apodization and the calibration process. 
Although we have included the imaginary parts of the data in our analysis we will discuss 
only the physically meaningful real results, while the imaginary results provide useful 
information on the errors. 
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4. Analysis Results 

4.1. The Galactic Plane 

In the Wright et al. (1991) analysis the dust continuum emission map of the sky 
was used as a spatial filter for averaging the emission from the spectral lines, effectively 
assuming that the spectral line intensities follow the continuum spatial distribution. That 
assumption is not made here. Rather, for each pixel we perform a simultaneous least 
squares fit of a continuum spectral "baseline" plus a series of spectral line profiles centered 
on the wavelengths of known lines. The result is a set of spectral line intensities for each 
sky position. We find that while the C + emission (Figure 1, top) follows the continuum 
closely, and the N + emission (Figure 1, bottom) follows it approximately, the C° and CO 
emission are distinct from the continuum template used by Wright et al., in that the C° 
and CO flux distributions are strongly concentrated towards the galactic center. 

In Figure 2a the longitude profiles of the C + and 158 fim far infrared continuum 
emission are shown for the entire galactic plane. The line and continuum emission are seen 
to track one another very well, as would be expected for PDR gas. The emission appears to 
peak towards the molecular ring, the galactic center direction, and in the spiral arms. Only 
towards the galactic center does the C + emission appear to be relatively weak compared 
with the far infrared continuum. In contrast, the 370 and 609 fim C° emission and the CO 
emission, shown in Figure 2 (i-p), peak strongly towards the galactic center direction. The 
J=2-l, 3-2, 4-3, and 5-4 lines are detected, while upper limits are placed on the CO J= 1-0 
and 6-5 transitions. 

Figure 2c and 2d show the profiles of the two N + transitions. Although the 122 /mi 
transition is noisier due to the reduced instrument sensitivity at that wavelength, the 
transitions clearly track one another, as expected. If we calculate a constant 122 to 205 fim 
intensity ratio we arrive at 0.9 ±0.1, consistent with the lower value given by Wright et al. 
(1991). However, the derivation of the intensity ratio of 0.9 is dominated by the stronger 
regions, and is not a particularly good fit overall. If, in addition to fitting for a constant 
ratio term, we include a term proportional to the 205 fim N + line intensity of the form 
7(N+ 122 //m)/7(N+ 205 fim) = 0.7 + 7(N+ 205 //m)/7 , then we find that the preferred 
solution is 7 = 7 X 10~ 5 erg s -1 cm -2 sr -1 . That is, in the regions where the 205 fim N + 
emission is strongest the intensity ratio is about unity, but in weaker regions the intensity 
ratio is indistinguishable from the 0.7 ratio expected for regions of low electron density. 

The spatial distributions of C + and N + are unmistakably related, but the N + intensity 
increases proportionally more than the continuum intensity while the C + intensity follows 
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the continuum in the plane of the Galaxy. Figure 3 shows that an excellent correlation 
exists between the fluxes of C + and N + , in the sense that /(N + ) oc /(C -1- ) 1 ' 5 , as pointed 
out by Petuchowski & Bennett (1993a, b) and Bennett & Hinshaw (1993). They suggest 
that the 1.5 power law may be the result of a volume to surface area ratio sampling of 
the interstellar gas. Figure 4 shows that J(C+) oc /(FIR) - 95 while J(N+) oc /(FIR) 1 - 53 , 
where the far IR intensity, /(FIR), is defined as the continuum intensity at each respective 
wavelength within the same effective bandwidth as the line. 

The 1.5 index power-law may be suggestive of a volume to surface area ratio geometric 
effect. In such a scenario the volumes of very large internally ionized regions give rise to 
N + emission. These regions are partially surrounded by neutral gas with the interface 
surface traced by C + emission. If the emission from each line of sight is dominated by 
one such region that is characterized by an effective radius R e ff } then /(N + ) oc /?e// an d 
/(C + ) oc Rljf, so from one line of sight to another the ratio will be characterized by 
/(N + )//(C + ) oc Riff- For the 1.5 power law to hold each line of sight must be better 
characterized by an effective radius than by the number of regions along each line of sight. 
If several regions contribute substantially to the emission along each line of sight, then 
/(N + ) would be linearly proportional to /(C + ). In this picture, the emission, which peaks 
in the direction of the molecular ring, is dominated by approximately beam-sized bubbles of 
ionized gas surrounded by neutral gas. For the COBE 7° beam, this corresponds to bubbles 
1 kpc in diameter. Comparable size bubbles are known to exist in external galaxies (e.g. 
Bruhweiler, Fitzurka & Gull 1991; Kamphius, Sancisi & van der Hulst 1991). Note that the 
N + emission need not arise from the volume of the bubble, but could be from swept up gas 
so long as the amount of gas is characterized by the bubble volume. 

The above scenario requires large-scale structures. Shibai et al. (1991) report results 
from a 3.4 arc-minute beam size survey of C + emission concluding that, "most [C II] 
emission of the Galactic disk is radiated not from discrete sources but from more extended 
and diffuse regions." Nakagawa et al. (1993) have shown C + maps made from observations 
with a 12 arc-minute beam that also shows structures that are much larger than the beam 
size. Doi et al. (1993) conclude that discrete sources are significant when looking through 
the length of the /=90° spiral arm. 

Massive and luminous and B stars tend to form in clusters and then become Type II 
supernovae that are correlated in space and time. The correlated supernovae create bubbles 
in the plane of the Galaxy. If a bubble is large enough it can break through the thickness 
of the disk. If the break-through is particularly energetic it will blow material into the 
halo of the Galaxy. These blowouts lead to Galactic fountains. In these models significant 
amounts of material from the galactic plane are thrown out of the plane and at least some 
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of this material rains back down onto the plane. These processes have been described by 
Bregman (1980), Houck & Bregman (1990), Mac Low & McCray (1988), Mac Low, McCray, 
& Norman (1989), Heiles (1990), Norman & Ikeuchi (1989), and Shapiro & Field (1976). 
Heiles (1990) estimates an area filling factor of breakthrough bubbles of 23% in our Galaxy. 
McKee & Williams (1993a, b) derive the somewhat lower filling factor of superbubbles of 
10%. We feel that the cycle involving giant molecular clouds, massive clustered and B 
stars, and clustered supernovae must be greatly enhanced in the molecular ring so that the 
filling factors of bubbles, quoted above, are probably underestimates for the ring. 

Petuchowski & Bennett (1993a, 1993b) examined the possible make-up of the COBE 
FIRAS line emission in terms of three morphologies: classical ionization bounded H II 
regions, boundary ionized clouds, and an extended low density warm ionized medium. They 
conclude that the classical ionization bounded H II regions contribute very little to the total 
line emission intensities, largely because of their small volume filling factors. The other two 
morphologies can occupy large filling factors and contribute significantly to the total signal 
intensity. Since bubbles, superbubbles, and the neutral material around their periphery, 
are examples of large filling factor regions, they are capable of dominating observed signals 
along lines of sight that also contain several, much smaller, classical H II regions. 

Volume-to-surface geometry is only one possible explanation for the 1.5 power 
scaling law between C + and N + . It may be that the C + scales proportionally with 
the far infrared continuum for the average PDR conditions in the Galaxy, but that 
J(N + ) oc /(FIR)^ 1 ' 5 oc /(UV)^ 1 ' 5 for some other physical reason, for example, having to do 
with the variations of color of the UV radiation field in the Galaxy. It is worth noting that 
a similar 1.5 power law relation has been reported for different transitions in an external 
galaxy (Petuchowski & Bennett 1993a, b). 

H I (c.f. Dickey & Lockman 1990) and C + appear to have different spatial distributions, 
but this does not contradict the picture that C + cools H I gas. Except where the emission is 
self-absorbed, the H I flux is directly proportional to the H I column density. On the other 
hand, the 158 fim C + flux does not, in general, scale directly with the gas density. Since 
C + emission is the result of a collisional process, one might naively expect C + emission to 
scale as the square of the gas density. In the low excitation PDR models of Hollenbach, 
Takahashi, & Tielens (1991), the 158 fim C + emergent intensity scales somewhat less than 
the square of the density for low densities and low UV field strengths. Nevertheless, since 
H I and C + intensities scale differently with density the fact that they have apparently 
different spatial distributions is not surprising. In addition, C + cools a fraction of the H 2 
gas (to a depth of Ay 4) as well as the H I gas, and warm H I gas has additional coolants 
available, further complicating comparisons of C + and H I distributions. 
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As discussed above, the 145.5 fim line of 0° is a potentially important PDR cooling 
line. Figure 2e shows the galactic plane distribution of the 145.5 fim line of neutral atomic 
oxygen, which Wright et al. (1991) listed as a nearly 2cr feature. The mean ratio to the C + 
intensity is J(C + )//(0°) = 71 ± 36, assuming a single intensity scaling factor. We caution 
that the 145.5 //m 0° line is close in wavelength to the very much stronger 158 fim C + line 
so that an error in our assumed instrument line response to C + at 158 fim of only 2% at 
145.5 fim could produce the 0° feature. 

Figure 2 also shows the galactic plane solutions for other selected transitions in the 
FIRAS wavelength range. 

4.2. The Galactic Center (< 500 pc) Region 

The physical properties of the Galactic Center (GC) region gas are unlike those 
elsewhere in the Galaxy. This is likely a result of the steep gravitational potential gradient 
within 500 pc of the GC where there are 10 8 Mq of molecular gas (Genzel 1992; Glisten 
1989). For molecular clouds near the GC to be stable against tidal disruption their mean H 2 
density must be about 10 4 cm -3 or more (Glisten & Downes 1980), in contrast to the typical 
disk molecular cloud density of 10 2 cm -3 (Scoville & Sanders 1987). Density-sensitive 
measurements of CS 1-0/2-1 show extensive high density regions within the central 500 pc 
(Tsuboi et al. 1989; Stark et al. 1989). The volume filling factor of molecular gas in the GC 
is two orders of magnitude higher than in the plane, resulting in mean gas densities of 100 
cm -3 in the GC versus 1 cm -3 in the plane (Sanders, Solomon, & Scoville 1984; Scoville & 
Sanders 1987). While individual clouds near the GC may be close to virialization, the larger 
cloud complexes are not (Stark & Bania 1986). The gas temperatures near the GC (40-150 
K) are much higher than the dust temperatures (< 30 K), as well as being higher than the 
temperatures of a typical disk molecular clouds (10 K) (Bally, Stark & Wilson 1987, Bally, 
Stark & Wilson 1988, Morris et al. 1983, Liszt & Burton 1978). Collisional heating of the 
gas by radiatively heated dust thus appears not to be the dominant heating mechanism 
near the GC. Blitz et al. (1985) report a decrease in the inferred cosmic ray rate toward the 
GC, based on COS-B 7-ray data, so a high rate of cosmic ray heating also appears to be 
ruled out. Heating by young stars appears to be ruled out as well because of the low dust 
temperatures and other indications of a low level of active star formation (Morris 1989; Cox 
& Laureijs 1989). The high gas densities themselves appear to rule out magnetic heating 
effects. The unusually large gas velocity-induced linewidths observed near the GC suggest 
the possibility of heating by turbulent energy dissipation (Wilson et al. 1982; Glisten et 
al. 1985). Since the turbulent energy dissipation time constant is short (~ 10 6 years), new 
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energy must be provided almost continuously, perhaps Irom the strong differential rotation 
in the central bulge (Fleck 1980). 

The AT&T Bell Labs CO J= 1-0 survey (Bally, Stark & Wilson 1987) shows a region 
that extends over galactic longitudes —1° < £ < +2° and about 0.4° of galactic latitude 
(i.e. about 400 pc X 50 pc in size). This region appears as a parallelogram in a £v diagram, 
indicating that much of the neutral material is in nearly rigid-disk rotation about the 
galactic center. Several H II regions detected by the H109a radio recombination line 
indicate star formation activity has occurred there in the last 10 8 years, and the positions of 
the H II regions suggest that star formation occurred in a 100 pc ring about the galactic 
center. 

Caution must be used in using CO to trace H 2 gas. For spherically symmetric 
neutral clouds in virial equilibrium with a radial velocity field (v oc r) the CO luminosity 
is given by Leo = (37rG'/4 / o) ' 5 TcoM (c.f. Scoville & Sanders 1987), where G is 
the gravitational constant, p is the cloud mass density, Too is the peak brightness 
temperature in the CO transition, M is the cloud mass, and Leo is in the usual units 
of K km s -1 pc 2 . The ratio of Leo to M is constant only as long as the geometry, 
dynamics, and composition of the clouds are constant, and Tco/p 1 ^ 2 is constant. Since 
temperature and density often rise together, this leads to a generally robust estimate that 
Xco = Nh 2 / Ico — 3.6 X 10 2 ° cm -2 (K km s -1 ) -1 . However, many assumptions go into the 
CO to H 2 conversion above. At one extreme, a molecular cloud can have all of its carbon in 
the form of C° and C + and no CO, but with H 2 at its core. In this case the CO intensity 
(zero!) underestimates the H 2 mass. At the other extreme a large dense molecular cloud 
with only a small surface thickness of C + can have almost all of its carbon in the form of 
CO. The fact that different volume-to-surface area ratios exist between molecular clouds in 
the disk and near the GC implies that one might expect a different CO to H 2 conversion 
factor. 

The far infrared fine structure lines, as seen by the COBE FIRAS, towards the GC 
are unlike any other region of the sky, and deserve special comment. 3.5° corresponds to 
500 pc at the GC, assuming a distance of 8.5 kpc. While there is a local maximum of 
emission of C + towards the GC, the C + emission has a global emission maximum towards 
the molecular ring. On the other hand, the CO and C° transitions peak strongly in the GC 
direction. Thus, the cooling of the GC neutral gas is being done proportionally more by the 
CO and C° forms of carbon than by C + . This would occur if the neutral gas towards the 
GC has a larger volume-to-surface area, or in other words, is better shielded, than neutral 
gas in the disk. Therefore, based on the observations of the COBE FIRAS carbon cooling 
lines, we suggest that there is an unusually large amount of well-shielded dense neutral gas 
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in the GC direction compared with other lines of sight. Since the GC neutral material is of 
higher density, relatively more cooling is expected from C° and CO than by C + because of 
the higher relative filling factor of high- Ay neutral gas. 

Alternately, differences in the FIRAS observations of the cooling lines towards the GC 
could be indicative of other effects such as a saturation of the 158 /tm C + line towards the 
GC or self-absorption of C + emission towards the galactic center, but the 1.5 power law 
scaling with the N + , discussed above, still holds in this region. 

The ground state Ooo — In 269 /tm transition of para-H 2 has never been detected 
astronomically. This is largely due to the enormous opacity of the Earth's atmosphere at 
this wavelength. The longitude profile of this water transition, shown in Figure 2h, displays 
an absorption-like dip towards the galactic center, although the signal-to-noise ratio is 
inadequate for us to be able to claim a detection of this line. The intensity dip at 269 /tm 
is ~ 2.5cr in frequency space. With about 5000 sky samples one expects 38 such dips on the 
sky at this wavelength and, in fact, there are 38. However, only towards the galactic center 
are three approximately independent dips so spatially close to one another. The error bars 
on the line intensity are larger near the galactic center because of uncertainty in the ability 
to calibrate and remove the larger continuum emission there. Thus we consider the water 
intensity dip suggestive, but not conclusive, and pose the question of whether such an 
absorption would be astrophysically reasonable. 

Water has a dipole moment of /t£> = 1.9 Debye, corresponding to the 269 /tm 
para- water transition rate A = 1.86 X 10~ 2 s -1 . If we assume the abundance of water is 
[H 2 0]/[H 2 ]= 10~ 6 , that 10% of all water resides in the gaseous ground state, and that the 
density of an absorbing cloud is typically 10 2 cm -3 , then the spectral line will be very 
optically thick. Hence, we expect optically thick water absorption in this transition to be 
common. For a linewidth corresponding to a velocity dispersion of 30 km s -1 (a typical 
width for lines near the GC), the FIRAS observation would suffer a frequency dilution 
factor of 150. Taking the mean value of the negative intensity dip as a real absorption 
gives A7 = 6 X 10~ 7 erg s -1 cm -2 sr -1 with a continuum intensity at that location and 
wavelength of I = 7 X 10~ 5 erg s -1 cm -2 sr -1 , so A7/7 = 9 X 10~ 3 . Correcting for the 
frequency filling factor gives 150A7/7 ~ 1, so the measurement is consistent with optically 
thick water absorption. The corresponding ground state transition of ortho-water is at 
538 /tm, also observable by FIRAS. This transition will also be optically thick, so any 
absorption feature in 269 /tm para-water should appear at 538 /tm at approximately a fifth 
of the 269 /tm intensity. There is no evidence for any absorption dip at 538 /tm (see Figure 
2g), but the expected level of A7 = 10~ 7 erg s -1 cm -2 sr -1 is below the level of the noise. 
Thus, while we can not claim a detection of this water line, the intensity dip observed at 
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269 /tm is consistent in wavelength, spatial position, and amplitude with what is expected 
astrophysically. 

The SubmiUimeter Wave Astronomy Satellite (SWAS) will observe the ortho- water 
ground state transition at 538 /tm, in addition to its analog in H^O at 547 /tm, 609 /tm C°, 
544 /tm 13 CO J=5-4, and 615 /tm O2. These lines fall in, or very near, the break between 
the high and low frequency FIRAS bands where the instrument response is greatly reduced. 
The 705.8 /tm transition of O2 is not detected by COBE FIRAS, and we place a limit on 
its intensity towards the galactic center of 2 X 10~ 7 erg s -1 cm -2 sr _1 . We place a limit of 
4 X 10~ 7 erg s -1 cm -2 sr -1 on the intensity of the 538.3 /tm transition of H 2 towards the 
galactic center. 

4.3. High Galactic Latitude Results 

Only the 158 /tm line of C + and the 205 /tm line of N + are sufficiently strong to allow 
useful estimates of their intensities at high galactic latitudes. 

We compute the ratio of the 158 /tm C + intensity to the H I column density for the 
entire sky at 7° resolution. Figures 5 and 6 show the relationship between the COBE 
158 /tm C + emission and H I column density (Dickey & Lockman 1990). A fit to the 
H I data gives a plane parallel approximation of N(H I) = 3.3 X IO 20 cm -2 csc|6| 
(Lockman, Hobbs, & Shull 1986). The difficulty in fitting for the C + cooling rate from 
the data in the low intensity portion of Figure 6 (top) is that not all of the C + emission 
comes from regions coextensive with atomic H I emission. C + emission also arises from 
molecular H 2 and ionized hydrogen regions. For regions where significant H 2 gas exists 
data points in Figure 6 are effectively shifted to the left, since the H I column density 
underestimates the total neutral hydrogen along the line of sight. For ionized hydrogen 
regions with significant C + emission, the data points in Figure 6 effectively take on 
a positive vertical displacement. Thus Figure 6 should contain an underlying linear 
relation with additional "contaminated" data above the linear relation. We can limit 
the contamination by selecting only data greater than 30° from the galactic center with 
N(HI) < 6 X 10 21 cm -2 . This eliminates most of the data where CO plays a major role in 
the cooling of the gas (i.e. where the hydrogen gas is in molecular H 2 form), and much of 
the gas where hydrogen is ionized in high column density lines of sight. A simple linear 
regression to the remaining data, shown in the lower panel of Figure 6, gives a cooling rate 
of 4:TrI(C + )/N(H I) = 2.8 X 10~ 26 erg s -1 per H-atom. However, these data remain clearly 
contaminated. Much of this contamination arises from the Orion and Gum Nebulae. We 
iteratively perform a median robust estimation of linear slope and intercept coefficients 
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and then remove the data that most deviates from the line. This results in a cooling rate 
of 2.65 X 10~ 26 erg s -1 per H-atom. The uncertainty in removing the "contamination" 
dominates, by far, the COBE FIRAS data uncertainties. We believe a "best" value for the 
diffuse C+ cooling rate is 4tt/(C + ) /N(H I) = (2.65 ±0.15) x 10" 26 erg s" 1 per H-atom. The 
contribution of C + emission from the ubiquitous Warm Ionized Medium can be estimated 
from the J(C + )-intercept of the plot in Figure 6. We place a limit of 



Bock et al. (1993) conclude that at a high galactic latitude the C + cooling rate is 
(2.6 ± 0.6) X 10~ 26 erg s -1 per H-atom in the small region of sky that they sampled and 
claim a detection of a nonzero intercept from ionized gas of 



The agreement of the cooling rates measured by Bock et al. with high resolution over a 
small sky region and our results with low resolution over nearly the full sky is especially 
reassuring. A UV-absorption derived C + cooling rate of 2.1 X 10~ 26 erg s -1 per H-atom 
was reported by Savage et al. (1993), less than the earlier result of Pottasch, Wesselius & 
van Duinen (1979) of 1.3 X 10~ 25 ergs s -1 per H-atom. Gry, Lequeux, & Boulanger (1992) 
reported a C + cooling rate of 3. 51^1 x 10~ 26 ergs s -1 per H-atom, but their suggestion that 
the COBE C + emission may come mainly from H II regions appears to be incorrect. 

The spin-orbit interaction removes the ground state degeneracy in singly ionized carbon 
causing the 158 fim magnetic dipole transition. Thus we consider C + as a two-level system 
and, denoting the 2 Pi/2 ground state as level "1" and the 2 P3/2 first excited state as "2", 
the critical density for level 2 is defined as n cr = A21/521, where q 2 \ = 8.629 X 10~ 6 £}/ g 2 T 1 / 2 
cm 3 s -1 is the collisional de-excitation rate and is a constant. For C + -e~ collisions we 
adopt = 0.444 (Hayes & Nussbaumer 1984), replacing the earlier laboratory result (1.34 
of Tambe 1977). For C+-H collisions we adopt = 0.00292 (Launay & Roueff 1977). The 
collisional excitation and de-excitation rates are related by q\ 2 = (#2/ ' 9\)<l2\e~ hv l kJ ■ The 
cooling function, Li (erg cm 3 s -1 ), resulting from collisions with partner "i" is given by 



U = {qi2)M = {g2/gi){q2i)Me- h ^ kT . 

If the 158 fim line originates in neutral gas then its intensity is given by (c.f. Crawford 
et al. 1985) 
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I(C + 158//m) = —hvA 
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where the summation runs over the various significant collision partners that excite C + 
(e.g. H I, H 2 , and electrons), hv = 1.26 X 10~ 14 ergs, the transition rate, A = 2.36 X 10~ 6 
s -1 , the statistical weights of the states are g = 2 J + 1 (g 2 = 4, §\ =2), and we assume 
a beam filling factor of = 1. iV#(C + ) is the column density of neutral hydrogen that 
is coextensive with the C + emission. The critical density, n cr depends on the the C + 
excitation partner. For collisions with H I, n cr (H I)~ 3000 cm -3 , and for collisions with H 2 , 
^cr(H2)~ 4000 cm -3 . For electron collisions n cr (e~)~ 10 cm -3 . If we assume that all of the 
emission results from coextensive H I and electron excitation in the Cold Neutral Medium 
(CNM), characterized by n = 30 cm -3 and T = 100 K, we can derive the C + cooling rate 
per H-atom. We can safely assume that n <C n cr in the CNM, so 



Using the definitions above relating the critical density with the cooling functions, we 
rewrite the C + cooling rate as 



where x e = n e /riH is the fractional electron density and the pressure, and p = nT is in the 
usual astrophysical units of cm -3 K. The temperature dependence of L/T is weak, making 
C + an excellent probe of the interstellar pressure in the CNM. The electron density results 
from both the photoionization of carbon and cosmic ray ionization of hydrogen, 



The intensity ratio is, I(C+ 158 /mi)//(iV + 205 /mi) = 19.6 ± 0.2, at |b| > 15°. Note that 
the predicted level for the CNM agrees well with the observations. The calculation can 
be repeated for the Warm Neutral Medium (WNM), characterized by n ~ 0.4 cm -3 and 
T = 7000 K. In this case the resulting intensity of C + emission is only a few percent of the 
CNM emission. This is not surprising since both the electron density and ion density are 
smaller, making collisional excitation much less effective. 

Finally, we consider the contribution of the Warm Ionized Medium (WIM) to 
the observed 158 fim C + intensity. This was computed by Reynolds (1992) to be 
Iwim{C + 158/xm) = 3.6 X 10~ 7 csc|6| erg cm -2 s -1 sr -1 . Bock et al. estimate that the 
WIM contribution to their particular sky region is (9 ± 3.5) X 10~ 8 erg s -1 cm -2 sr -1 . We 
derive an upper limit using an extrapolation to N(H I)=0 of the fit to the data in Figure 




(10) 
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6, Iwim(C + 158/rni) < 6 x 10 7 erg s 1 cm 2 sr 1 (95% CL). We conclude that the high 
galactic latitude C + emission predominantly arises from the CNM. 

We have examined a grid of model parameters to estimate the properties of the high 
latitude gas cooled by the C + emission. There are many ranges of the parameters that fit 
the data well. For example, T = 100 K, n = 10 cm" 3 , ^ = 3.8 x 10" 4 , and ( CR = 1 x 10" 17 

s -1 corresponds to a pressure of 1000 cm -3 K. T = 100 K, n = 20 cm -3 , = 2.0 X 10~ 4 , 
and (cr = 5 X 10~ 18 s -1 corresponds to a pressure of 2000 cm -3 K. T = 80 K, n = 15 cm -3 , 
= 2.7 X 10~ 4 , and (cr = 1 x 10~ 16 s -1 corresponds to a pressure of 1200 cm -3 K. 

When Nh > 6.5 X 10 21 cm -2 the slope in Figure 6 dramatically increases. Assuming the 
C + cooling rate remains at (2.65 ± 0.15) X 10~ 26 erg s -1 Hatom -1 , and using the estimated 
total C + luminosity from our Galaxy from COBE of 10 7 ' 7 solar luminosities (Wright et al. 
1991), we derive that the 158 fim line from C + is cooling > 10 9 solar masses of gas. 

In Figure 7 we plot the C + line-to-continuum ratio as a function of galactic latitude. 
This ratio is defined using the instrument resolution bandwidth for the continuum flux. 
The line-to-continuum intensity ratio falls from about unity on the plane, to about one 
half at |&| ~ 50°, to the point where C + is undetectable (< 10% of the plane value) to the 
COBE FIRAS instrument above \b\ ~ 65°. We interpret these data as indicating that the 
volume filling fraction of excited surfaces of cold neutral gas that preferentially give rise to 
C + emission decreases with increasing galactic latitude. If this is the case then the high 
galactic latitude far infrared emission that remains must be coming almost entirely from 
warmer gas. The higher excitation conditions in the warmer gas allow many atomic levels 
to become populated, and thus several higher energy cooling transitions can dominate over 
the ground state transition of C + . Also, the warm gas is much less dense than the colder 
medium, and thus the efficiency of 158 fim cooling is greatly reduced. It is difficult to 
determine what fraction of high latitude far infrared emission comes from warm neutral 
versus warm ionized gas. This question was addressed by Boulanger & Perault (1988). 
They correlated H I emission with a 100 fim IRAS map and found a very high degree of 
correlation, and concluded that the far IR is coming from neutral regions. However, they 
also point out that their result leaves no room for the expected IR emission from the Warm 
Ionized Medium. This can be reconciled, however, if the H I and ionized gas are at least 
partially correlated with one another. Results of a survey using the Wisconsin Ha Mapper 
(WHAM), now under development, will help to resolve these questions. 



4.4. Total Ionizing Photon Rate in the Galaxy 
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The total ionizing photon rate in the Galaxy, </>, can be estimated from the total 
luminosity of the 205 //m N + observations from the COBE FIRAS, since the Galaxy is 
transparent at this wavelength. We revise the estimated luminosity of the 205 //m N + 
observations from Wright et al. (1991) upward by 30% to account for the fact that that 
the N + intensity follows the 1.5 rather than 1.0 power of the continuum, as was assumed in 
that work, now giving -Zv 2 o5,umNii = 2.6 X 10 40 erg s -1 inside of the solar circle. The ionizing 
photon rate is given by 

^ _ -^205 AtmNII [H + ] OiH+ 

/L/m Nil [N + ] C^205 (im Nil 

where hu 205limNII = 9.7 x 10" 15 erg, [N+]/[H+] = 6.8 x 10" 5 , and ^os^mNii = 
6.7 x 10- 8 (T/8000 K)-°- 5 n e cm 3 s" 1 and a H + = 6 x 10- 13 (T/8000 K)"°- 5 n e cm 3 s" 1 are 
the recombination rates. The result is <f> = 3.5 X 10 53 ionizing photons per second with 
~50% uncertainty. This is in remarkable agreement with Mezger (1978), who arrived at 
3.0 X 10 53 s -1 with 20-30% uncertainty. Our derived rate can be thought of in terms of the 
ionization provided by 25,000 O-stars throughout the Galaxy, each with an ionizing photon 
rate of 1.4 X 10 49 s" 1 . 



5. Conclusions 

1. The COBE FIRAS has conducted an unbiased survey of the far-infrared emission from 
our Galaxy. The first results of this survey were reported by Wright et al. (1991). We 
report here the morphological distribution along the galactic plane (b = 0°) of the spectral 
line emission for all of the detected spectral lines. 

2. We find that the strongest spectral line is the 158 fim fine structure transition from C + . 
This emission peaks in the molecular ring and can be seen to follow the spiral arm structure 
of the Galaxy. The intensity of the C + line scales linearly with the far-infrared continuum 
in the plane of the Galaxy. 

3. The 205 //m line of N + generally follows the 158 fim C + line distributions, but the 
intensities scale as I(N + 205 /im) oc I(C + 158 /mi) 1 ' 5 in the inner galactic plane. It has 
been suggested that the 1.5 power law may indicate a volume to surface area morphological 
sampling. 

4. We present the cosecant dependence of the 158 fim fine structure transition from C + as 

I(C + 158 /mi) = (1.43 ± 0.12) x 10 -6 esc |b| erg cm" 2 s" 1 sr" 1 . 

We estimate the C + cooling rate to be (2.65 ±0.15) X 10~ 26 erg s -1 Hatom -1 , and conclude 
that the emission is arising almost entirely from the Cold Neutral Medium with a pressure 
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of about 1000-2000 cm -3 K. We find that the C + emission falls off much faster than the far 
IR continuum with increasing galactic latitude, indicating that the high galactic latitude 
far IR arises from a warm medium. 

5. We present the cosecant dependence of the 205 fim fine structure transition from N + as 

J(N + 205 /mi) = (4 ± 1) x 10" 8 esc |b| erg cm -2 s _1 sr _1 

and conclude that the emission is arising entirely from the Warm Ionized Medium. 

6. We estimate the total ionizing photon rate in the inner Galaxy to be <f> = 3.5 X 10 53 
ionizing photons per second, with ~50% uncertainty, based on the 205 fim N + transition. 
This is in remarkable agreement with Mezger (1978), who arrived at 3.0 X 10 53 s -1 with 
20-30% uncertainty. 

We are grateful to Eli Dwek, Mike Hauser, Carl Heiles, Chris McKee, Bill Reach, Ron 
Reynolds, and Sam Petuchowski for useful discussions. COBE is supported by NASA's 
Astrophysics Division. The Goddard Space Flight Center (GSFC), under the scientific 
guidance of the COBE Science Working Group, is responsible for the development and 
operation of COBE. We are grateful for the contributions of the entire COBE Science Team 
and the support personnel at Goddard's Cosmology Data Analysis Center (CDAC). 
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FIGURE CAPTIONS 

Figure 1: (Color Plate) The maps are projections of the full sky in galactic coordinates. 
The plane of the Milky Way is horizontal in the middle of the map with the galactic center 
at the center. Galactic longitude i = 90° is left of center. The maps are smoothed to 
10° resolution, (top): A full sky map of 158 /tm C + emission from the COBE FIRAS 
experiment, (bottom): A full sky map of 205 /tm N + emission from the COBE FIRAS 
experiment. 

Figure 2: Longitude profiles of the intensity solution for spectral lines at given 
wavelengths for the entire galactic plane. The intensity in each bin represents an average 
over 5° in galactic longitude and ±5° in galactic latitude. Error bars include detector noise 
and gain uncertainties, which affect the accuracy of the continuum removal. There are very 
few measurements in the £ = 90°, 6 = 0° bin, so we caution against overinterpretation of 
results based only on this data point, (a) The detected C + 2p: 2 P 3 / 2 — 2 P\/2 fine structure 
ground state transition at 158 /tm. The dashed line is the far infrared continuum emission 
at that wavelength and within the same effective bandwidth, (b) The undetected O2 3 2 — I2 
at 705.8 //in. (c) The detected N + 2p 2 : 3 Pi — 3 P fine structure ground state transition 
at 205.3 /im. The dashed line is 0.07 times the far infrared continuum emission at that 
wavelength and within the same effective bandwidth, (d) The detected N + 2p 2 : 3 P 2 — 3 P\ 
transition at 121.9 /tm. (e) The possibly detected (see text) 0° 2p 4 : 3 P — 3 Pi transition 
at 145.5 //in. (f) The undetected Si 3p 2 : 3 P 1 - 3 P Si transition at 129.7 /mi. (g) The 
undetected ortho-H 2 lio — loi ground state transition at 538.3 /tm. (h) The possibly 
detected (see text) para-H 2 In — loo ground state transition at 269.3 /tm. (i) The 
undetected CO J= 1-0 transition, (j) The detected CO J=2-l transition, (k) The detected 
CO J=3-2 transition. (1) The detected CO J=4-3 transition, (m) The detected CO J=5-4 
transition, (n) The undetected CO J=6-5 transition, (o) The detected C° 2p 2 : 3 Pi — 3 P 
transition at 609.1 /tm. (p) The detected C° 2p 2 : 3 P 2 — 3 Pi transition at 370.4 /tm. 

Figure 3: An excellent correlation exists between the intensity of the C + and N + 
intensities in the Galactic plane, in the sense that 7(N + ) oc /(C -1- ) 1 ' 5 . 

Figure 4: J(C+) oc I(FIR) - 95 while J(N+) oc I(FIR) 1 - 53 in the Galactic plane. 

Figure 5: (Color Plate) The maps are projections of the full sky in galactic 
coordinates. The plane of the Milky Way is horizontal in the middle of the map with 
the galactic center at the center. Galactic longitude £ = 90° is left of center. The maps 
are smoothed to 10° resolution, (top): The full sky distribution of neutral hydrogen H I 
emission, described by Dickey & Lockman (1990). (bottom): The C + cooling rate per 
hydrogen atom, 47r/(C + 158 fj,m)/N(H P). 
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Figure 6: (top): The observed 158 //m C + intensity versus the H I column density, 
N(H) cm -2 , as described by Dickey & Lockman (1990). (bottom): Data restricted to at least 
30° Irom the galactic center and N(H) < 6 X 10 21 H-atoms cm -2 . The three lines correspond 
to our best estimate of the C + cooling rate of (2.65 ± 0.15) X 10~ 26 erg s -1 Hatom -1 and 
its 95% confidence upper and lower bounds. This assumes that the H I measurements are 
accurate and that C + and H I emission are coextensive. 

Figure 7: The C + line-to-continuum intensity ratio as a function of galactic latitude. 
The far IR continuum is defined as the best fit continuum intensity within the same effective 
bandwidth (i.e. instrument resolution bandwidth) as the C + line emission. The C + emission 
decreases with increasing galactic latitude until it becomes undetectable at \b\ > 65°. 
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